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Abstract 

We report studies of the decay — > DcpK^, where Dqp denotes neutral D mesons that decay 
to CP eigenstates. The analysis is based on a 29.1 fb~^ data sample of collected at the T(4S') 
resonance with the Belle detector at the KEKB asymmetric e+e" storage ring. Ratios of branching 
fractions of Cabibbo-supprcssed to Cabibbo- favored processes involving Dcp are determined to be 
B{B- DiK-)/B{B- DiTT-) = 0.125±0.036±0.010 and ^(5" ^ D2K-)/B{B' D2Tt-) = 
0.119 zb 0.028 lb 0.006, where indices 1 and 2 represent the CP = +1 and CP = —1 eigenstates of 
the — D^ system, respectively. We also extract the partial rate asymmetries for B^ — DcpK^, 
finding = 0.29 ± 0.26 ± 0.05 and A2 = -0.22 ± 0.24 ± 0.04. 

PACS numbers: 12.15.Hh, 13.25.Hw 
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Direct CP- violation in DcpK^ decay, where Dcp denotes neutral D mesons that 

decay to CP eigenstates, provides a promising way to extract the angle 03 of the Cabibbo- 
Kobayashi-Maskawa unitarity triangle [|I], 0]. A partial rate asymmetry between the DcpK~ 
and DcpK^ final states can arise from interference between b ^ c and h ^ u processes 
shown in Fig. |^. The relation of the partial rate asymmetry Acp and the 03 angle |^ is 
given by: 



B{B- Di^2K~) + B{B+ Di^2K+) 

2r sin 6' sin 03 
1 + + 2r cos 6' cos < 



(1) 



^3 



where indices 1 and 2 denote the CP = +1 and CP = —1 eigenstates of the neutral D 
mesons, r is the ratio of the amplitudes, r = 1^4(5^ D^K^)/A{B~ D^K~)\, and 5' is 
the strong phase difference between the two amplitudes, with 5' = 5 for Di and 5' = 5 + tt 
for D2- This asymmetry can have a non-zero value when both 03 and 5 are non-zero. In 
principle, one can constrain the angle 03 from the measurement of asymmetries Ai and A2- 

Experimentally, B DK processes have been studied using measurements of the ratio of 
the Cabibbo-suppressed process B~ D^K~ to the Cabibbo-favored process B~ D^tt~ . 
Belle [§ measured R = B{B- D^K-)/B{B- D^n-) = 0.079 ± 0.009 ± 0.006, while 
CLEO 1^ reported R = 0.055±0.014±0.005. Assuming factorization, R is naively expected 
to be tan^ Oc{fK/fw)^ ~ 0.074 in a tree-level approximation, where 9c is the Cabibbo angle, 
and fx and /^r are the decay constants. Both measurements are in agreement with this 
theoretical expectation. 

In this Letter, we report the first measurement of B^ — > DcpK^ decay. We also give 
a comparison of the ratio of branching fractions R for the fiavor specific state and the 
CP = ±1 eigenstates, and a determination of the asymmetries Ai^2- The results are base on 
a 29.1 fb~^ data sample collected on the T(4S') resonance with the Belle detector [Q at the 
KEKB asymmetric e+e~ collider 0. This corresponds to approximately 31.3 million BB 
events. The inclusion of charged conjugate states is implied throughout this Letter unless 
otherwise stated. 

Belle is a general-purpose detector with a 1.5 T superconducting solenoid magnet that can 
distinguish the Cabibbo-suppressed process B~ D^K~ from the Cabibbo-favored process 
B~ — ^ D^TT~ by means of particle identification and kinematic separation. A charged-particle 
tracking system, covering approximately 90% of the solid angle in the center-of-mass (cm) 
frame, is comprised of a three-layer silicon vertex detector (SVD) and a 50-layer central 
drift chamber (CDC). Identification of charged hadron species is accomplished by combining 
responses from silica aerogel Cerenkov counters (ACC), the time-of-fiight detector (TOP) 
and the dE/dx measurement from the CDC; it provides more than 2. Sex K/n separation 
for laboratory momenta up to 3.5 GeV/c. A CsI(Tl) Electromagnetic Calorimeter (ECL) 
located inside the solenoid coil is used for photon (7) detection. A detailed description of 
the Belle detector can be found elsewhere 0. 

We analyze both B~ D^K~ and B~ D^n^, which are collectively referred to as 
B^ D^h'^ decays. The processes B^ D^K~ and B^ D^n^ are distinguished by a 
tight requirement on the particle identification of the prompt {K~ or vr") and the effect 
of the mass difference at the final stage of the event selection. The decay B~ — > D^tt~ is used 
as a control sample to establish constraints on kinematic variables, resolution of detectors, 
evaluation of systematic uncertainties and normalization of results. 
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Flavor specific meson candidates are reconstructed via Df ^ K~ti^\ for CP = +1 
eigenstates we use Di — > K^K^ and 7r~7r"'" and for CP = —1 we use D2 Ksn^, Ksuj, 
Ks(p, KsTj and KsTj' . 

Tlie Ks TT+TT^ candidates are reconstructed from two oppositely charged tracks with an 
invariant mass within ±3(j of the nominal Ks mass. Candidate 7r° mesons are reconstructed 
from pairs of 7's, each with energy greater than 30 MeV, and are required to have a recon- 
structed mass within ±3ct of the nominal vr" mass. We form candidate rj and 77' mesons using 
the 77 and rjn^TT' decay modes with mass ranges of 0.495 GeV/c^ < M(77) < 0.578 GeV/c^ 
and 0.904 GeV/c^ < M(r77r+7r^) < 1.003 GeV/c^, respectively. The r] momentum is required 
to be greater than 500 MeV/c. 

Candidate uj and mesons are reconstructed from tt+tt^tt'^ and K^K~ combinations 
with invariant masses in the ranges 0.733 GeV/c^ < M(7r+7r^7r'^) < 0.819 GeV/c^ and 
1.007 GeV/c^ < M{K^K~) < 1.031 GeV/c^, respectively. A helicity angle requirement for 
the UJ and (p mesons of | cos 6hei\ > 0.4 reduces the contributions from non-resonant — >■ 
/TsTT+Tr^TT^ and KsK^K~ decays to a negligible level. For Ksto, the — > K*~ 

background is explicitly rejected by vetoing any Kst^~ combination that forms an invariant 
mass within ±75 MeV/c^ of the nominal K* mass. 

For each charged track from the meson decay, the particle identification system is 
used to determine a K/vr likelihood ratio P{K/-n) = Ck/{^k + ^n), where Ck and are 
kaon and pion likelihoods 0]. For K~tt^ , K^K~^, and 7r~7r"'", the charged tracks with 

P{K/ti) < 0.7 are assigned to be pions, while kaons are required to satisfy P{K/tt) > 0.3 
for D° ^ K-7r+ and P{K/n) > 0.7 for D° R- . Pions from candidate D° tt-h+ 
decays are required to have momentum greater than 0.8 GeV/c, and the candidate is vetoed 
if either pion, when combined with any other track in the event, has an invariant mass that 
is within ±50 MeV/c^ of the nominal J /if) mass, or within ±20 MeV/c^ of the nominal 
mass. 

Candidate mesons are also required to have an invariant mass within ±2.5cr of 
the nominal mass, where a is the measured mass resolution, which ranges from 5 to 
18 MeV/c^, depending on the decay channel. Tracks and photons from the candidate 
final states, except for Kstt^, KgTj and Kstj', are refitted according to the nominal me- 
son mass hypothesis and the reconstructed vertex position, in order to improve momentum 
determination. 

We reconstruct —>■ Dh~ events using the laboratory constrained mass (M/c). Mi^ is 
the B candidate mass calculated from the laboratory momenta by using an e^e~ BB 



hypothesis: Mi^ = \J [E^^'^Y ~ (Pb)"^, where ps is the B candidate's laboratory momentum 
vector and E''£^ = ± P^e ■ Pb), where s is square of the cm energy, is the 

laboratory momentum vector of the B meson candidate, and Pee and E^e are the laboratory 
momentum and energy of the e~^e~ system, respectively. Mic is independent of the mass 
assumption of the particle to be used, and is suitable for both B~ Dti~ and B~ —>■ DK~ 
simultaneously. We accept B candidates with 5.27 GeV/c^ < Mic < 5.29 GeV/c^. 

Background events from e+e^ — > qq continuum processes are rejected using event shape 
variables that distinguish between spherical BB events and jet-like continuum events. We 
construct a Fisher discriminant, = Y.i=2,a(^iRi° ± J2t=i , where ai, Pi are optimized 
coefficients to maximize the discrimination between BB events and continuum events, and 



Rf", Rf° are modified Fox- Wolfram moments [|r^ |T2|. Furthermore, the variable cos 63 is 



used, where 6b is the angle between the B fiight direction in the T(4S) rest frame and the 
beam axis. A single likelihood variable is formed from the probability density functions of 
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T and cos 9b. The likelihood ratio is defined as CTZ = £sig/('Csig + '^cont) where £sig and 
£cont are likelihoods defined by and cos 9b for the signal and continuum background, 
respectively. Since each sub-decay mode has different backgrounds, we optimize the cut 
of the likelihood ratio for each mode by maximizing S/ ^/S + N, where S and N denote the 



numbers of signal and background events estimated by a Monte Carlo simulation |10|. For 
the —>■ K~7i~^ mode, the requirement CTZ > 0.4 retains 87.1% of signal and 26.4% of the 
continuum background. 

To distinguish between D^K^ and D^n^ processes, we form two samples by 

requiring that the prompt hadron h~ satisfies either P[K/'n-) > 0.8 {B^ D^K~ enriched 
sample) or P{K/'k) < 0.8 {B~ D^7t~ enriched sample). Finally, the AE distribution, 
calculated by assigning the pion mass to the prompt hadron h~ , is examined to distinguish 
between B~ — > D^K~ and B~ —>■ D^tt" processes, where AE is the energy difference in 
the cm frame: AE = EJ^ + E"^ - E^^m- The signals for B' D'^n' and ^ D^R- 
peak in the AE distributions at MeV and —49 MeV respectively, as can be seen in Fig. |^. 
In the B~ —>■ D^K~ enriched sample, we can also see a peak around AE = MeV due to 
misidentified pions from B' — > D^tt~ in Fig. @(b), (d), and (f). 

The numbers of B~ D^tt~ and B^ — D^K~ events are extracted by fitting dou- 
ble Gaussian functions with different central values and widths to the AE distribution. 
Backgrounds originate from qq continuum events and BB events. Continuum events 
are distributed over the entire AE region, and the shape of this background is deter- 
mined by fitting a linear function to the AE distribution in the Mi^. sideband region 
(5.20 GeV/c^ < Mi^ < 5.26 GeV/c^). BB backgrounds, such as ^ D°p" and B D*n- 
processes, are mostly seen at negative values of AE; MC simulations are used to obtain the 
shape of their distribution. 

In the fits to the B^ — > D'^tt^ enriched AE distribution, the signal peak position and 
width, and the normalization of continuum and BB backgrounds are free parameters. On the 
other hand, we calibrated the shape parameters of the B^ D^K~ signal using the B^ — »• 
D^Ti~ data following the procedure described in |Q, which accounts for the kinematical shifts 
and smearing of the AE peaks caused by the incorrect mass assignments of prompt hadrons. 
For the B~ —>■ D^K~ fit, the shape parameters of the double Gaussian of the signal are fixed 
at values determined from the fit to the B^ D^tt^ sample. The peak position and width of 
the B D^K~ signal events are determined by fitting the B^ D^7i~ distribution using 
a kaon mass hypothesis for the prompt pion, where the relative peak position is reversed 
with respect to the origin. The shape parameters for the feed-across from B~ D°7r~ is 
fixed by the fit results of the B~ — ^ D^vr" enriched sample. 

The fit results are shown as solid curves in Fig. 0. The number of events in the D^K~ 
signal and D^7r~ feed-across, and the statistical significance of the D^K~ signal are given in 



Table |. The statistical significance is defined as 2 ln(£o/'Cmaa:), where Cmax is the max- 
imum likelihood in the AE fit and Cq is the likelihood when the signal yield is constrained 
to be zero. The statistical significance of both DiK~ and D2K~ signals is over 5. Oct. 
Experimentally, the ratio of branching fractions is determined as follows: 



N{B- D^K-) r]{B- D 



TT e TT 



N{B-~^DH-) 7]{B-^DH<-) e{Ky ^' 

where N is the number of events obtained, r] is the signal detection efficiency, and e is the 
prompt pion or kaon identification efficiency. The signal detection efficiencies are determined 
from MC simulation: ti{B' D^K~) is approximately 5% lower than ri{B~ —>■ D^7f~) 
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due to kaon decays in flight. Kaon and pion identification efficiencies e{K) and e{n) are 
experimentally determined from a kinematically selected sample of D*^ D^ir^ and 
— > K~7r^ decays with K~ and tt"*" mesons from candidates that are in the same 
cm momentum (p'^™) and polar angle regions as prompt hadrons in the — > Dh~ decay 
(2.1 GeV/c < p^™ < 2.5 GeV/c). With our requirements, P{K/'k) > 0.8 gives a kaon identi- 
fication efficiency of e{K) = 0.778 ±0.005 with a pion misidentification rate of 0.024 ±0.002, 
and P{K/tt) < 0.8 gives a pion identification efficiency of e(7r) = 0.972 ± 0.007. 

Since the kinematics of the B~ D^K~ and B~ — > D'^tt' processes are similar, the 
systematic uncertainties from detection efficiencies cancel in the ratios of branching fractions. 
The dominant systematic errors are the uncertainties in the shapes of backgrounds in the 
AE distributions (5.1% — 7.9%), which are determined by varying the background shape 
of the fitting function by ±la, and K/ir identification efficiencies (1.2%). The sum of the 
uncertainties for each mode are combined in quadrature to determine the total systematic 
errors for the ratios. 

The resulting measurements of R are listed with their statistical and systematic errors 
in Table |. These are the first observations of the decays B~ —>■ DcpK~. As a check, the 
result for the flavor specific decay B~ — > D^h^ , K~tt^ is listed as well, and is found 

to be consistent with previous measurements. We find no deviation of the R ratio for the 
B^ DcpK~ processes from the corresponding flavor specific modes beyond statistical 
errors. 

Direct CP violation is investigated by measuring the partial rate asymmetries in B"^ — >■ 
Di 2K'^ decays. We fit the yields of B~^ and B~ events separately for each mode and deter- 
mine the partial rate asymmetries Ai^2 are shown in Table The partial rate asymmetry 
is consistent with zero for the flavor specific mode, which is expected to have no asymmetry. 
To construct 90% confidence intervals of Ai and A2, we combine statistical and systematic 
errors in quadrature, and assume that the total error is distributed as a Gaussian. We find 
—0.14 < ^1 < 0.73 and —0.62 < A2 < 0.18, consistent with zero asymmetry. 

The main sources of systematic errors for the partial rate asymmetries Ai^2 are asym- 
metries in the measured background (1.5% — 3.9%), intrinsic asymmetry in the Belle detec- 
tor (3.6%), and kaon identification (1.0%). We observe 1154.6 ± 35.4 B+ D°7r+,Z)° ^ 
K^n^ candidates and 1073.5±34.5 B^ D^tt' , K~tt^ candidates. This is consistent 

with our expectation that the detector has no significant intrinsic charge asymmetry. Using 
MC simulation, we find that the contribution of the non-resonant contaminations (~ 0.1%) 
of u and can be neglected. 

In conclusion, using 29.1 fb~^ of data collected with the Belle detector, we report studies 
of the decays B^ DcpK^, where Dcp are the neutral D meson CP eigenstates. This 
is the first stage of a program to measure the angle 03 in the CKM unitarity triangle. 
The ratios of branching fractions R for the decays B~ DcpK~ and B~ DcpT^~ are 
consistent with that for the flavor specific decay within errors. The measured partial rate 
asymmetries Ai^2 are consistent with zero within large errors. In the future, with 300 fb~^ of 
data, an accuracy of better than 0.1 in the measurements of i?i_2 and Ai^2 in DcpK~ 
is expected and we will begin to constrain the angle 03 . 

We wish to thank the KEKB accelerator group for the excellent operation of the KEKB 
accelerator. We acknowledge support from the Ministry of Education, Culture, Sports, Sci- 
ence, and Technology of Japan and the Japan Society for the Promotion of Science; the Aus- 
tralian Research Council and the Australian Department of Industry, Science and Resources; 
the Department of Science and Technology of India; the BK21 program of the Ministry of 



7 



Education of Korea and the CHEP SRC program of the Korea Science and Engineering 
Foundation; the Pohsh State Committee for Scientific Research under contract No.2P03B 
17017; the Ministry of Science and Technology of Russian Federation; the National Science 
Council and the Ministry of Education of Taiwan; the Japan- Taiwan Cooperative Program 
of the Interchange Association; National Science Foundation of China under Contract No. 
10175071; and the U.S. Department of Energy. 



* on leave from Nova Gorica Polytechnic, Nova Gorica 
[1] M. Gronau and D. Wyler, Phys. Lett. B265 172 (1991). 

[2] I. Dunietz, Phys. Lett. B270 75 (1991); D. Atwood, I. Dunietz and A. Soni, Phys. Rev. Lett. 

78 3257 (1997). 
[3] H. Quinn and A.I. Sanda, Eur. Phys. J. C15 626 (2000). 
[4] K. Abe et al. (Belle Collaboration), Phys. Rev. Lett. 87 111801 (2001). 
[5] M. Athanas et al. (CLEO Collaboration), Phys. Rev. Lett. 80 5493 (1998). 
[6] A. Abashian et al. (Belle Collaboration), Nucl. Instr. and Meth. A479, 117 (2002). 
[7] KEKB B-factory Design Report, KEK Report 95-7 (1995), unpublished. 
[8] The charge conjugated mode(-D'^ — > K^n^) is implicitly included. 

[9] The uj helicity angle is defined as the angle between the flight direction in the T(4S)'s 
frame and the normal to the uj decay plane in the w's rest frame. The 4> helicity angle is defined 
as the angle between the (p flight direction in the D's rest frame and the direction of the (/>'s 
daughter in the i;^'s rest frame. 

[10] We use the QQ i?-meson decay event generator developed by the CLEO Collaboration 
( Pittp:/ /www. Ins. cor nell.edu /public/CLEO/soft/QQ) and GEANT3 for the detector simu- 
lation; CERN Program Library Long Writeup W5013, CERN, 1993 

[11] K. Abe et al. (Belle Collaboration), Phys. Rev. Lett. 87, 101801 (2001). 

[12] The index so(oo) indicates that the moment is calculated from pairs of particles where only 
one (neither) particle comes from the B candidate. 




8 




AE(GeV) AE(GeV) 



FIG. 2: IS.E distributions for B D^tt /K candidates and fit results: (a) B Dfir , 
(b) B- DfK-, (c) B- Diir-, (d) B^ DiR-, (e) ^ Davr-, and (f) B' D2K~, 
where in each case the pion mass is assigned to the prompt ■k~/K^. Dotted (dashed) lines show 
the distributions of DK{D-jt) signals. The shaded plot shows the continuum background and the 
remaining component from BB background is estimated and fitted by Monte Carlo simulation. In 
the DK plots, the dashed curves show the D'^tt feed-across. 
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TABLE I: Results of fits for the D 7r~ and D K~ decay modes. The event yields, the feed-across 
from £)°7r~ to the D^K~ signal region, statistical significance of D^K~ signals, efficiencies and 
branching fraction ratios (i?) are given. Efficiencies are determined by weighting according to the 
measured sub-components (77 = Z^j??iH(D° — Xi)). 

B- D^K- B- r»°7r- Stat. B' D^tt' Efficiency(%) branching fraction ratio 

events feed-across sig. events ri{D^7r~)/r]{D^ K~) R 

B- Dfh- 161.7 ± 14.5 51.3 ± 9.7 16.9 2245.1 ±51.0 1.703/1.639 0.094 ± 0.009 ± 0.007 

B- ^ Dih- 22.9 ±6.1 9.6 ± 4.4 5.1 240.1 ± 16.7 0.173/0.165 0.125 ± 0.036 ± 0.010 

B- ^ D2h- 26.1 ±6.5 4.9 ±4.1 5.5 290.6 ± 19.1 0.184/0.173 0.119 ± 0.028 ± 0.006 



TABLE IL Summary of measured partial rate asymmetries. 





Mode 


N{B+) 


NiB-) 


Acp 


90% C.L. 




DfK^ 


80.6 ± 10.1 


81.1 ± 10.4 


0.003 ± 0.089 ± 0.037 


-0.15 <Af< 0.16 




DiK^ 


8.1 ±3.9 


14.7 ±4.6 


0.29 ± 0.26 ± 0.05 


-0.14 <Ai< 0.73 




^D2K^ 


16.4 ±5.5 


10.6 ±4.2 


-0.22 ± 0.24 ± 0.04 


-0.62 <A2< 0.18 
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